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Introduction: Understanding variations in activity 
with composition is an essential step in improving pre- 
diction of partition coefficients during magma evolu- 
tion. Variations in activity with composition are com- 
plex and do not generally exhibit ideal behavior rela- 
tive to a traditional melt-component set. Although 
deviations from component ideality can be modeled 
numerically by simply fitting to compositional vari- 
ables (such as in a regular or subregular solution 
model), such models have not been particularly suc- 
cessful for describing variations in trace component 
activities. A better approach might be to try to identify 
components that do a better job of describing the be- 
havior of the species in the melt. 

A significant part of the nonideal contribution to 
activity is due not to intrinsic properties of the compo- 
nent in the melt, but to our ignorance about the state 
and mixing properties of the component in the melt. In 
general, the more realistically that a component’s en- 
ergy state and mixing properties are portrayed, the 
smaller the nonideal component becomes. Looked at 
another way, this means that the variation in activity 
coefficient with composition provides insight into the 
mixing properties of the component in the melt. 

We have measured the activity of R 112 O 3 in a vari- 
ety of melt compositions and propose that variations 
arise due to contributions from both a Ru 3+ and a 0 2 - in 
the melt. The first of these two contributors, a Ru 3 +, can 
be understood as depending on the mixing between 
Ru and A1 , Mg , and Fe . Variations in the sec- 
ond contributor, a 0 2 -, can be modeled by a buffering 
reaction with bridging and non-bridging oxygens. 



Si0 2 

AI203 

MgO 

CaO 

run T, C 

diop 

56.0 

0.1 

17.8 

26.7 


gg 

48.4 

8.0 

25.7 

16.1 

1580 

gg2%AI-B 

55.2 

2.4 

29.9 

11.6 

1578 

gg4%AI-B 

53.5 

4.6 

29.5 

11.2 

1587 

gglOCa-B* 

53.3 

9.6 

30.5 

6.2 

1587 

gg14Ca-B* 

51.8 

9.3 

29.6 

8.9 

1586 

gg17Mg 

55.8 

9.4 

16.5 

18.6 

1580 

gg22Mg 

51.4 

8.7 

21.2 

17.8 

1584 

gg46Si 

45.5 

8.3 

27.3 

17.7 

1580 

gg52Si* 

51.7 

7.6 

25.4 

15.0 

1580 


Tablet: Experimental compositions (wt%). 0.3wt% RU 2 O 3 


Electrochemical Measurement of Ru 2 0 3 activi- 
ties: Electrodes were inserted into silicate melt beads 
of various compositions (Table 1) suspended on Pt- 
wire loops in a 1-atm gas mixing furnace. An electri- 
cal potential was imposed between the electrodes, the 
imposed potential increasing along a step ramp with a 
pulse imposed on each step (Fig. 1). Current flows 
between electrodes when electroactive species in the 
melt are oxidized or reduced at the electrodes. The 
resulting current was measured at the top and bottom 
of the voltage pulse, and the difference (the differential 
current) was plotted against potential. The peak of the 
resulting curve is related to the activity coefficient for 
the particular electroactive species (Ru 2 0 3 ) in the melt 
[1,2,3]. 



Fig. 1 Potential step function for electrochemistry experiments. 
ph=pulse height (ph), 

The Ruthenium peak occurred on the shoulder of a 
large interfering peak associated either with oxidation 
of silicate melt, or oxidation of the Pt electrode; the 
latter being the most likely. This large peak was not 
very reproducible, and therefore subtracting it out as 
background to yield a net Ruthenium peak was not 
successful. 

Instead, we took advantage of the observation from 
theoretical curves that when two different curves are 
added together, a sharp slope change occurs near the 
potential at which the two curves contribute equal 
amounts to the total current (Fig. 2). By averaging the 
E value at which this slope break occurs with the E 
value at V 2 that current, we derived a value for the Ru 
peak (Fig. 2). The resultant peak had an uncertainty of 
about ±14mV. 

An activity coefficient (relative to a standard state 
of diopsidic melt) can be calculated from the expres- 
sion: 

lnymeit = -nF/RT(Epk diop - Epk melt ) + ln[(molar vol- 
ume diopsidic melt)/(molar volume melt of interest)]. 
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(where Epk is measured at the same f0 2 in each com- 
position, n=number of electrons in the reaction). 
Results are given in Table 2. 


gg 

2.05 

gg2%AI-B 

3.02 

gg4%AI-B 

2.82 

gglOCa-B 

3.38 

gg17Mg 

1.97 

gg22Mg 

2.33 

gg46Si 

2.46 

gg52Si 

4.20 

Diop 

1.00 


Table 2: Measured values for yRu 2 0 3 relative to 
a standard state of melt of diopsidic composition. 

Modeling activity for Ru0 l s and the oxide ion: 

Activities of oxides in silicate melts can be understood 
as depending on the separate variations in activities of 
the cation and the oxide ion [3], after a relationship of 
the sort: 

RuOi.s = Ru 3+ + 3/20 2 ' 
and 

yRuO l5 = Yr u 3+ ■ (o 2 )) 32 

The activity of the oxide ion, O 2 ", can be consid- 
ered to be buffered by reactions of the type: 

Si-O-Si + O 2- — > 2(Si-0)" 

Therefore, the activity of the oxide ion can be de- 
fined by 

a 0 2 " = K*(Si-0) 2 /(Si-0-Si) 

If we interpret Si-O to be non-bridging oxygens 
(NBO), and Si-O-Si to be bridging oxygens (BO), and 
define a G 2 ' to equal 1 in diopsidic melt, then 
a 0 2 " = 0.75 • (NBO) 2 /(BO) 

where NBO = [2*0 -4*(Si0 2 +0.5*A10i. 5 )]/0; BO 
= 1-NBO. 

Inferring the mixing behavior of Ru ,+ in silicate 
melts from statistical correlations: A value for y Ru 3+ 
can be calculated from the expression: y Ru 3+ = 
yRu0i. 5 /(0 2 ") 3/2 . Insofar as RuOi.s behaves ideally 
relative to this proposed component set, a Ru 3+ cc 
X Ru 3+/(Z cations with which Ru’ mixes randomly) and 
y Ru 3+ oc 1/(1 cations with which Ru 31 mixes randomly). 
Therefore, a regression to the expression 
l/y Ru 3+ = A + B • AIO 1.5 + C • MgO + D • CaO 
provides an indication of the relative degree with 
which Ru 3 1 mixes with different cations in the melt. 

The regression to 9 observations yields the follow- 
ing: R 2 = 0.76; A=-4.57(1.09); B=4.58(2.20); 



Fig. 2 Two slopes are drawn on the composite curve as 
shown. This break in slope occurs near the point where 
the two curves contribute equal amounts to the composite 
curve. The peak of the smaller curve (related to Ru 2 0 3 ) is 
found by averaging E at the slope break with E at 'A the 
current at the slope break 

C=9.27(2.03); D=12.7(2.5). 

For comparison purposes, an analogous regression 
for electrochemically derived FeO activity data (8 ob- 
servations) to the expression 

l/y Fe 2 + = A + B • A10 15 + C • MgO + D • CaO 
yields the results 

R 2 =0.76 A= -1.67(1.41); B=0.43(2.95); C=4.64(2.11); 
D=5. 91(3.56). 

The primary difference between the Ru 31 and Fe 2+ 
results is that the A 1 0 1 . 5 coefficient for Ru 1 is larger 
than for Fe 2+ (relative to the other coefficients) and is 
more statistically significant. This is consistent with 
the trivalent Ru 3+ being more likely to substitute for 
Al 3+ on melt sites than is the divalent Fe 2+ . This sug- 
gests that simply assuming that network-modifiers 
substitute on network modifying sites and network- 
formers on network forming sites is inadequate. 

Conclusions: Ru’ 1 mixes on melt sites containing 
Al’ + , Mg 2+ , and Ca 2+ . It’s mixing behavior differs 
from Fe 2+ , indicating that a simple model in which 
network-modifying cations mix equally on network 
modifying sites is not correct. Rather, different cations 
mix differently on melt sites according to their size and 
electrical charge. Correlations between activity coeffi- 
cients and composition can provide insight into the 
melt components with which Ru’ 1 (or other cations) 
mix. 
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